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Abstract 



We show that the eigenvalues of the first order partial differential equa- 
tion derived by quasi-classical approximation of the Schrodinger equation 
can be computed from the trace of a classical operator. The derived trace 
formula is different from the Gutzwiller trace formula. The Fredholm de- 
terminant of the new operator is an entire function of the complex energy 
plane in contrast to the semi-classical spectral determinant derived from 
the Gutzwiller trace formula. 

1 Introduction 

The Gutzwiller trace formulajlj is the most important tool for semiclassical 
quantization of non-integrable systems. The trace formula, in general, does not 
converge close to the semiclassical energy eigenvalues or resonances where we 
would like to use it. The spectral determinant, derived by Vorosj| from the 
trace formula, has much better mathematical properties. However it also does 
not converge on the whole complex energy plane || and cannot be used to de- 
termine all the complex resonances of a system. One possible solution of the 
convergence problem is to introduce smeared trace formulas or spectral determi- 
nants. These methods limit the accuracy of the resonances and allows only qual- 
itative comparison of the exact quantum and the semiclassical spectra. Trace 
formulas and spectral determinants can be introduced for classical systems H 




including thermodynamic formalism B and their analytical properties can be 
studied. The convergence problems of these formulas and the semiclassical con- 
vergence problems have the same origin H . For a large class of classical systems 
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with 'nice' mathematical properties the recent theorem of H. H. Rugh|| assures 
the convergence of classical spectral (Fredholm) determinants for a large class 
of classical systems with 'nice' mathematical properties. This theory can be ex- 
tended to semiclassical trace formulas ||, and the semiclassical resonances can 
be computed without convergence problems. In addition, the numerical error of 
the semiclassical resonances decreases super-exponentially ~ exp(— n 13 ) {(3 > 1) 
as a function of the maximal topological length n of periodic orbits. 

The theory in Ref.|| has been developed on an abstract mathematical 
ground. Here, the connection of this mathematical description and the quasi- 
classical approximation, introduced by Maslovjl^], is discussed. The term quasi- 
classical is more appropriate than semiclassical since the Maslov type descrip- 
tion leads to a pure classical Perron-Frobenius operator in a natural way. The 
Fredholm determinant of this operator can serve as a well defined and mathe- 
matically clean starting point of the semiclassical approximation of the spectra. 
One additional advantage of this description is that the wave function evolves 
along one single classical trajectory and we do not have to compute sums over 
increasing numbers of classical trajectories as in computations involving Van 
Vleck formula |0J. 

2 The quasi-classical approximation 

Following mostly Ref. , we give a summary of the quasi-classical approxima- 
tion, which was worked out by Maslov|l(| in this form. This approximation was 
used earlier by Liouville, Green, Stokes, Rayleigh and others. 

The Schrodinger equation for a single particle of unit mass in a potential U 

is 

ih^ = ~^ + U{q)^ (1) 

where A is the Laplace operator and ip(q,t) is the wave function. The Ti is the 
small parameter of the problem, and the quasi-classical approximation gives the 
solution of this equation for h — > 0. This asymptotic solution is of short wave 
form 

i>{q,t) = v(q,t)e lS ^/ h , (2) 
where the amplitude Lp(q,t) and the phase S(q,t) are smooth real functions on 
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some bounded region of the configuration space. Substituting (j^) in the Schro- 
dinger equation, we can derive partial differential equations (PDE) for the am- 
plitude and the phase with initial conditions ip(q, 0) = <po{q) and S(q, 0) = So(q) 
respectively. In the limit h — > omission of the non-leading terms proportional 
to h 2 leads to the following quasi-classical partial differential equations 

^ + div( Q VS) = 0, (3) 

— + H(q,\7S) = 0, (4) 

where g(q,t) — ip 2 (q,t) and H(q,p) is the Hamilton function. Equation (||), 
which is the Hamilton-Jacobi equation, is a first order partial differential equa- 
tion whose solution corresponding to an initial S(q,0) = So(q), can be given 
independently from (|J). However, the continuity equation (||) can be solved 
only if the solution S(q, t) of (^) is known. In other words, the Hamilton-Jacobi 
equation is autonomous while the continuity equation is driven by the solution 
S(q,t) of @. 

We know from the theory of first order PDE's[jl2| that their solutions lead to 
ordinary differential equations (ODE). As is well known, the Hamilton-Jacobi 
equation leads to the Hamilton differential equations 

. = dH(q,p) 
dp 

where the new variable 

p = VS(q,t) 

has been introduced. In the PDE description we evolve the whole function 
S(q, t) and compute its gradient at a given point go- Computation of the gradient 
requires information about the behavior of the function S(q, t) in the vicinity 
of go and can not be recovered from the value of S(qo,t) alone. In the ODE 
description we evolve both q and p = VS(q, t). Therefore at a given time, we do 
not have to compute the gradient from S(q, t) since the evolution is local in the 
(q,p) space. From the ODE description we can reconstruct the whole solution 
of the PDE as 

S(q', t) = S(q 7 0) + / L(q(T),q(r))dT, (6) 
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where we have to integrate the Lagrange function along the phase space trajec- 
tory with 



q' = q(t),q = q(0),VS(q', t) = p(t),VS(q, 0) = p(0). 

For a generic S(q, 0) only one such trajectory exists. 

The 'local' solution of the continuity equation (|J) is also straightforward. 
At a given starting point go the momentum is given by po — VSo(go), an d 
the amplitude of the wave function is fo(q)- After time t the coordinate go 
evolves to q(t) and po to p(t) according to the full Hamiltonian flow. The 
new amplitude can be derived from the probability conservation as follows: 
Take an infinitesimal initial d dimensional directed volume V(qg) around go in 
the configuration space. The points of this volume evolve to the infinitesimal 
directed volume V(q(t)) around q(t) according to the Hamiltonian flow. The 
value of the amplitude changes according to 

where the sign ±1 reflects the ambiguity of the transformation from the density 
g to if. Careful analysis shows that the minus sign refers to the case where the 
orientation of the final volume is the opposite the initial one. The ratio of vol- 
umes is independent of the way we specify the initial infinitesimal volume. To 
understand which quantities determine the volume ratio, we specify an initial 
directed parallelepiped around go with edges given by d independent infinites- 
imal vectors <5qi, <5q2, . . . <5qd- After time t these vectors are transformed into 
(5q' 1; <5q' 2 , ■ ■ • Sq 1 d by the flow. The initial point (go,Po) and the set of initial vec- 
tors do not specify uniquely the image vectors. This is a consequence of the fact 
that the infinitesimal vectors are transformed by the full Jacobian matrix of the 
Hamilton flow, and we have not specified the infinitesimal momentum vectors 
around po corresponding to these vectors. The initial function So(q) determines 
a set <5pi, <5p2, ■ ■ ■ Spd of vectors through the second derivative matrix: 

5p = MSq, M, i = ^^, (8) 

which we shall call the curvature matrix. The initial curvature matrix M is an 
important initial condition, and we are not able to compute the image of the 
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volume without it. The vector (5qj,#pj) is transformed by the Jacobi matrix 

Sq't = JqqSqi + JqpSpi, (9) 

Sp\ = JpqSqi+JppSpi, (10) 

(11) 

where the subscripts q and p denote the corresponding [d x d] submatrices of 
the full [2d x 2d] Jacobian. The Jacobian is determined by the initial condition 
(qo,Po) and can be computed as a time ordered integral along the phase space 
trajectory 

J( 9j p,i) = Texp{jf dTT> 2 H(q(r),p(T))} , (12) 

where D 2 H(q,p) denotes the second symplectic derivative matrix of the Hamil- 
tonian and T is the time ordering operator. The curvature matrix of the function 
S(q, t) around q(t) defines a (^|) type relation between the infinitesimal vectors 

*p' = M'*q', ^ = ^P . (13) 
dq J dq l 

Using (||) and ([TJ), we can eliminate the vectors Spi and Sp^ and can get 
relations between the initial and final infinitesimal configuration vectors and 
the curvature matrices 

6q! = (3 qq +3 qp M)Sci, (14) 
M' = (J^+Jp^CJ^+J^M)" 1 . (15) 

From the first relation (|TJ) we can read off the volume ratio: 

V(a') 

-^l=det(J qq +J qp M). (16) 

The second relation (|l^) is a recursion relation for M and can be considered as 
the matrix generalization of the usual fractional rational transformation. From 
( [l5| ) we can derive a differential equation for M(t), if we substitute the elements 
of the infinitesimal Jacobi matrix. This differential equation 

m= tit + M 7fir + 7nr M + m tit m - 17 

\Oqoq opoq oqop opop J 

is a driven one since the second partial derivatives of the Hamilton function 
should be computed along the phase space trajectory. If we solve this differential 
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equation along the phase space trajectory the volume ratio can be written as a 
time ordered integral along the phase space and M(i) trajectory 



V(q ) 



= Texp 



Tr 



d 2 H d 2 H 



dpdq dpdp 



M 



(It 



The square root of the volume ratio is also a time ordered integral: 



( v(q') 

\V(qo) 



-1/2 



T exp 



Tr 



d 2 H d 2 H • 
dpdq dpdp 



(18) 



(19) 



The computation of this expression requires some care when the solution of the 
differential equation ( |l~7| ) is singular. Close to a singularity, where 

M(t -> t c ) = oo, 



(20) 



we can neglect the non-leading terms from fll7|) and use the solution of 



M = M 



d 2 H 
dpdp 



M. 



The second derivative matrix can be decomposed into combinations of dyads 
and their eigenvalues in the usual way 



d 2 H 
dpdp 



^AiD i; D,D, ,),,T>, 



(21) 



The solution close to the singularity can be a linear combination of some of 
these dyads corresponding to singular directions I: 

R 

(22) 



R 

M(t) = V^D ; , 



where R is the number of singular directions. The time ordered integral close 
to the singularity is dominated by 



( v(q(t c +o)) 
U(s(*io)) 



-1/2 



= exp 



te+O 



R 



tr 



-rfr 



(23) 



This integral can be computed by adding infinitesimal imaginary value ie to t 
and taking the e — > limit 

fV{q{t%,))^- 1/2 



\V(q(t c _ )) 



exp(w(fl/2)). 



(24) 
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Between two singular points the time ordered integral is positive and gives the 
absolute value of the volume ratio. Notice that R counts the number of rank 
reductions of the matrix M along the classical path, and it is also a function on 
the initial condition Mo. 

Now we have everything needed to describe the time evolution of a quasi- 
classical wave function. The wave function at time t is now 



where go is the starting point of a classical trajectory with initial momentum 
V<So(<Zo) which ends up in q after time t with momentum \7S(q,t), and the 
volume ratio is determined by the curvature matrix M = didjS(q,0). 

3 Time evolution a la Maslov 

At this point it is possible to express the volume ratio and the momentum with 
the second and first derivatives of the minimal action between q' and q. In 
this way we recover the usual Van Vleck propagator Q. The spirit of equation 
( p5| ) is that the wave amplitude ip at time t and at coordinate q' is determined 
by the amplitude at t = at coordinate q. In calculations involving the Van 
Vleck operator kernel this nice property is lost, and we have to compute lots of 
trajectories to compute the volume ratio and we have to know the whole initial 
wave function too. However we have a better option. We can keep track of the 
variables p and M along only one trajectory and compute (^|) and the volume 
ratio (|l|). This means that the evolution takes place on the extended (q,p, M) 
space. We can introduce classical density functions ^ defined on this space. The 
wave function then corresponds to the special function 



The evolution of a general classical density function on the extended space 
according to (|||) can be rewritten in terms of a classical transfer operator 



^(g'.j/.M'.t) = / dqdpdMC(q',p',M',t\ q,p, M, 0)ip(q,p, M, 0), (27) 



i } (q' 7 t) = ^q' 7 t)e^'^ h = ± 




Mlo)e tSoM/h (25) 




(26) 
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with the kernel 

e l ™+Ja dT ^ + ^ TT {^ + ^ M } 6(q' -q\q,p))5( P ' -p\q,p))5(M' -M\q,p,M)), 

(28) 

where q t (q,p), p t (q,p) and M t (q,p 7 M) denote the evolution of q, p and M 
from the initial coordinates q,p = \7So(q) and M = didjSo(q) during the time 
t, and v = N + R/2. The time ordered integrals should be computed along 
the full trajectory, and also the number of rank reductions R and the number 
of orientation changes N. The Dirac deltas assure that the operator connects 
coordinates, which are connected by the classical dynamics, and give the correct 
amplitude. This operator can evolve densities, which are not of the form (p6|), 
and therefore we can expect that only a part of its spectrum has relevance to 
semiclassics, but all the semiclassical eigenvalues will be contained in its spectra. 



4 Entire Fredholm determinants 



The spectral or Fredholm determinant of the operator (28) can be defined by 
the identity 

det(l -C) = exp (- ^ ^Tr£"^J , (29) 

where the traces of the powers of the operator are 

Tr£" = J dqdpdMC n (q,p,M\ q,p,M). (30) 

The zeroes of the Fredholm determinant on the complex energy plane yield the 
eigenvalues of the operator. The simplest application of ( p8| ) is to 3-dimensional 
hyperbolic Hamiltonian flows. In this case the spectral or Fredholm determinant 
of the operator (E8) is given byM 



I 1 e irS p (E)+iv p ir 

Mi-*) - H-E r | A ;|(i-i/Ap 2 A - 



I A r 11/2 |A r h 5 / 2 

Ap ' r " 1 - I/A* + 1 - 1/Af ' [ ' 

where p is the index of primitive periodic orbits, A p is the corresponding eigen- 
value of the stability matrix, and r is the repetition number. The practical 
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advantage of (|T]) over the more familiar Gutzwiller-Voros and Ruelle type zeta 
functions was demonstrated by detailed numerical studies (?J of the related quan- 
tum Fredholm determinant |Q . It can be shown |l5| that the Fredholm deter- 
minant obtained by keeping only one of the terms in the sum in (|^) is an entire 
function in the whole energy plane. This enables us to show that the Gutzwiller- 
Voros spectral determinant for Axiom A flows is meromorphic in the complex 
E plane, as it can be written as a ratio of entire functions^. The non-physical 
eigenvalues of (28) can be removed. 



5 Conclusions 

In conclusion, we have constructed a classical evolution operator for semi- 
classical quantization based on Maslov's quasi-classical quantization method, 
and derived a new determinant for periodic orbit quantization of chaotic dy- 
namical systems. The main virtue of the determinant is that the theorem 
of H.H. RughQ, applicable to multiplicative evolution operators such as (p8|), 
implies that this determinant should be entire for Axiom A flows, i.e. free of 
poles in the entire complex energy plane. Our numerical tests of the three disk 
system also support the above claims flEf . 
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